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Abstract. The integrated mean curvature of a simplicial manifold is well understood 
in both Regge Calculus and Discrete Differential Geometry. However, a well motivated 
pointwise definition of curvature requires a careful choice of volume over which to 
uniformly distribute the local integrated curvature. We show that hybrid cells formed 
using both the simplicial lattice and its circumcentric dual emerge as a remarkably 
natural structure for the distribution of this local integrated curvature. These hybrid 
cells form a complete tessellation of the simplicial manifold, contain a geometric 
orthonormal basis, and are also shown to give a pointwise mean curvature with a 
natural interpretation as a fractional rate of change of the normal vector. 
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1. Introduction 

In representing the curvature of a manifold in some ambient space, the extrinsic 
curvature gives a measure for the deviation of normal vectors to the manifold about 
a given point. This deviation can be measured in a number of ways, including: (?) the 
difference between a parallel transported normal from the normal at a nearby point, 
(??) the fractional rate of change of the area of a small region along a one-parameter 
family of manifolds orthogonal to the normal vectors, or (???) the Lie derivative of the 
intrinsic metric along the normal, for the same one-parameter family. Each of these 
three definitions yields the same extrinsic curvature tensor on taking the infinitesimal 
limit for a smooth manifold. 

In both Regge Calculus (RC) and Discrete Differential Geometry (DDG), smooth 
manifolds are approximated by piecewise-linear simplicial manifolds. These are formed 
by connecting blocks of flat space, with each block given by an n-simplex. For each of 
these, the interior geometry is Euclidean, or Minkowski for 4-simplices in RC, which 
was developed specifically to give a discrete version of General Relativity (GR). These 
n-simplices are then connected by identifying their (n — l)-dimensional faces. 

Llnlike the continuum case, the different ways of defining curvature on a piecewise- 
flat manifold do not necessarily agree. However, the integral of both the intrinsic and 
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extrinsic curvatures, over certain paths, tend to be well-defined and consistent with 
the continuum. In the original work of Regge [T], it was shown that the integrated 
Gaussian curvature around co-dimension 2 simplices is all that is required to give an 
action corresponding with the Einstein-Hilbert action of continuum GR. Any volume 
that is used to distribute the integrated curvature over, to give a field tensor, cancels 
and so the action is invariant to such a choice. A similar integrated mean extrinsic 
curvature was first found by Hartle & Sorkin [2], in extending the Regge action to 
include boundary terms. These integrated curvatures have also been studied in the 
DDG community, see (3j and references therein, with the expression for the integrated 
mean curvature appearing as far back as 1840 in the work of Steiner [4j. 

Though the integrated curvatures are all that is required for defining the action in a 
vacuum spacetime, a pointwise curvature is necessary for the inclusion of matter. Such 
a definition requires the choice of some set of volumes to distribute the local integrated 
curvature over. This leaves room for ambiguities to arise. When the simplicial manifold 
forms a Delaunay triangulation, it’s dual Voronoi lattice is formed by joining the 
circumcenters of n-simplices. This circumcentric dual has long been known as providing 
one of the most natural dual areas for the intrinsic curvature [HUH! [7J IE]. More recent 
work has used hybrid volumes, depending on both the Delaunay and Voronoi lattices, 
to define pointwise scalar curvature and pointwise Riemann and Ricci curvature tensors 
PHD]. This has even led to a discrete version of the Ricci flow mi. 

Here it is shown that a particular type of hybrid volume, constructed using 
both the simplicial lattice and its circumcentric dual lattice, gives a remarkably rich 
interpretation for a pointwise extrinsic curvature. The geometry of these volumes 
is completely determined by the edge-lengths of the simplicial manifold, they give a 
complete tessellation of the manifold and each contains a natural geometric orthonormal 
frame. The distributed mean curvature can also be viewed as a fractional change in 
normal over a specific displacement. 

In what follows, the theory is first developed for the most simple case of a 2- 
dimensional surface in E 3 . The details of a simplicial submanifold are outlined, and an 
expression for the extrinsic curvature associated with each edge is derived in section [2] 
The circumcentric dual and hybrid volumes are then introduced in section P along with 
the new distributed curvature and its fractional interpretation. More general extensions 
are considered in section [4] with the intrinsic properties of simplicial manifolds of higher 
dimension described in section 14.11 defining the circumcentric dual and hybrid volumes. 
Hypermanifolds of E n are specifically treated in section 14.21 with the interpretations 
from the 2-dimensional case requiring little alteration. Generalizations to arbitrary 
dimensional submanifolds, and submanifolds of non-Euclidean manifolds are then brief 
outlined in section l4~3l 
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2. Simplicial Surfaces in 3-Dimensional Euclidean Space 

A piecewise-flat simplicial surface S, embedded in a 3-dimensional Euclidean space M, 
is a surface which can be decomposed into a collection of flat triangles (2-simplices). 
The flat geometry of each triangle is completely determined by the length of its edges 
l, which are given by the Euclidean 3-metric. This allows the internal geometry of the 
surface to be completely determined by the lengths of all of the edges. 

Since the boundary edge between two adjacent triangles has the same length in 
each triangle, and the geometry internal to each triangle is Euclidean, the geometry 
from one triangle can be smoothly developed across an edge to any adjacent triangle. 
Intrinsic deviations from Euclidean 2-space only occur when a closed path contains one 
or more vertices. The intrinsic curvature is then constructed from the deficit angle e 
around a vertex, given by the difference between the sum of the angles at the vertex 
and a full rotation in Euclidean space, see figure [1] 




Figure 1 . A simplicial surface consisting of 5 equilateral triangles (a), and the 
development of one of the triangles (b), giving the deficit angle e around the central 
vertex. 


On a smooth surface in E 3 , the extrinsic curvature can be defined as the 
infinitesimal parallel transport of the unit normal h, projected onto the surface. The sign 
convention of Misner, Thorne & Wheeler da will be used here. As shown in figure 
integrating the curvature along a path on the surface gives the finite parallel transport 
of hi, projected onto the surface at the point x 2 : 


2(a) 



I<; dx b 


sin 6v a ~ 9v a . 


( 1 ) 


The unit vector v a lies on the surface, in the direction of the projection of the vector 
hi — h 2 at the point x 2 . Notice that if the surface in figure 2(a) is concavely curved, 
the vector v a will point in the opposite direction. The small angle approximation for 
sin 6 can be taken whenever the curvature changes slowly with respect to the distance 
between the two points. 

On a piecewise-flat simplicial surface, the normal vector remains unchanged across 
a given triangle. The extrinsic curvature is thus concentrated on the edges l, or ‘hinges’, 
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between any pair of adjacent triangles. This is similar to the concentration of the 
intrinsic curvature on the vertices. To find the extrinsic curvature associated with an 
edge l, an arbitrary convex region D is defined, enclosing l and entirely contained within 


the two triangles separated by /, see figure 2(b). Within D, an orthonormal coordinate 
basis (x,y) is also defined, with y parallel to the edge l. For a given value of y, the 
integral of along a path 7 (y), parallel to the x coordinate from one boundary to 
the other, is given by the difference 9 between the normal vectors to each of the two 
triangles: 


K°: dx = 0v a = Ox 


h(y) 


( 2 ) 


The variation between the two normals can only occur perpendicular to the ‘hinge’ l, 
so the unit vector v a must be parallel to the x-coordinate. The sign is then given by 
the orientation of the projection of hi — relative to the x-direct ion. This will be 


positive for a concave, and negative for a convex hinging (as in figure 2(b)). The small 


angle approximation here simply assumes that the lattice is dense enough to give a small 
variation in the normals between adjacent triangles. 




(a) Smooth Integrated Curvature 


(b) Simplicial Integrated Curvature 


Figure 2. Integrated curvature along a finite path 7 on a smooth surface (a), and 
across an edge l of a simplicial surface (b). 

Since the y coordinate is parallel to /, integration of the extrinsic curvature along y 
leaves the normal fixed, and so the y components of can be assumed to vanish. The 
mean curvature K, given by the trace of the curvature tensor, is therefore equivalent to 
the xx component. For a dense enough lattice, the integral of K along x is then given 
simply by the angle 6 between the normal vectors of the two triangles, for any y: 

[ K dx = f Kf dx = 0, (3) 

J'V(y) •''Y (y) 

with 9 again positive for a concave curvature, and negative for convex (as in the figures). 
The integrated mean curvature over the entire region D can now be found, since the 
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integral of K along the x coordinate is invariant to the choice of y: 


rl r 


KdA = 


'D 


K dx 


h(y) 


dy = 9 dy = l 6 . 


( 4 ) 


Notice here that the label l is used simultaneously to represent both a particular edge 
and the measure of that edge. This convention will be used for a number of labels, 
though it should be unambiguous which use is intended in each situation. 

The value of the integrated mean curvature 19 does not depend on the shape of the 
region D , as long as it is convex, contains all of l, and does not cross any other edge. It 
is therefore possible to reduce the size of the region on either side of Z, with the result 
remaining unchanged in the limit as D converges to the edge l itself. In this case, the 
integrated mean curvature associated with an edge can be defined as: 


IMCi := 19 , (5) 

depending simply on the length of the edge, and the angle between the normals of the 
two triangles on either side. This expression was derived from the boundary of the 
action by Hartle & Sorkin [2], and by Brewin ra using a similar argument to here, but 
employing a family of smooth surfaces limiting to the simplicial surface. In the context 
of Discrete Differential Geometry, this expression can also be found in [3], where it was 
credited to Jakob Steiner from 1840 [1]. 


3. Circumcentric Dual and Mean Curvature Density 

To find a distribution of the curvature tensor over the entire surface, the circumcentric 
dual lattice leads to a remarkably natural definition for the region D[. This lattice 
consists of lines within the simplicial surface, joining the circumcenters of adjacent 
triangles, as illustrated in figure [3] For two triangles separated by an edge Z, the line 
joining the circumcenters of the two triangles can be shown to be perpendicular to l, 
intersecting it at its midpoint. This circumcentric line is considered dual to the simplicial 
edge Z, and denoted A/. Since the circumcenters of each triangle can be determined from 
the edge lengths of the triangle, the length of a dual edge A/ is determined entirely from 
the five edges of the two triangles it connects. 

A hybrid cell Hi is defined as the quadrilateral with diagonals given by Z and A i. 
These cells form a complete tessellation of the simplicial surface. The internal geometry 
of each hybrid cell is Euclidean, formed by the development from one of the triangles 
intersecting the cell, across the edge Z. One of the biggest advantages of the circumcentric 
hybrid cells, comes from the orthogonality of Z and A;. Unit vectors in these directions 
form a natural orthonormal basis for Hi, defining a piecewise-linear Cartan frame on S. 
The area of a hybrid cell is also conveniently given by the equation: 

Hi = — l A; , 

t 2 

due to the orthogonality of Z with its circumcentric dual A/. 


( 6 ) 
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(b) Circumcentric Hybrid Area 


Figure 3. The development of the simplicial surface from figure |T| with the 
circumcentric dual edges shown (a), and the hybrid cell formed by the edge l and 
its circumcentric dual A i (b). 


Defining the arbitrary region Di by the circumcentric hybrid volume Hi leads 
naturally to a distributed mean curvature, spread evenly over Hg. 

IMQ l e 


K, = 


H, 


l A, 


e 


( 7 ) 


The average distance from one side of Hi to the other, orthogonal to l , is given by 
half the length of the dual edge A/. This gives an alternative view of this distributed 
curvature, as the change in the normal vector from one side of Hi to the other, divided 
by the average distance across the cell: 


0 _ 9 2 9 

dave 2 ^ 

This gives an analogue of the continuous notion of the curvature as the rate of change of 
the normal vector with respect to a displacement along the surface. The role played by 
this average distance is again a special case of the circumcentric hybrid cells, resulting 
from its linear boundaries, and A i as the maximum extension of Hi orthogonal to the 
edge l. 

The mean curvature associated with a triangle itself can now be found. A first 
requirement is that the integrated mean curvature over all of the triangles of S, is 
equivalent to the sum of the integrated mean curvatures associated to each of the edges: 


IMC = Y. K ‘ A ‘ = Y. IMC ‘ • ( 9 ) 

t l 

with the subscripts t representing each triangle in the simplicial manifold. Using the 
densitised mean curvature given by the circumcentric hybrid cells: 




V imc, 

l 

E K ‘ H ‘ 
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where dt represents the boundary of the triangle t, and Hi\ t represents the restriction 
of the hybrid cell Hi to that part which lies in the triangle t. These new cells Hi\ t are 
known as reduced hybrid cells. The sum over all triangles t can now be dropped from 
both sides, with the integrated mean curvature for a given triangle given by: 

IMC t = K t A t = Y, K > H ’\t , (11) 

ledt 

seen as an average of the distributed curvatures Ki for the hybrid cells intersecting the 
triangle t, weighted according the areas Hi\ t of the restriction of their associated hybrid 
cells to t. With an expression for Ki from (JTJ) and OH]), and since Hig is simply a triangle 
itself: 

iMc t = ■ (12) 

ledt ledt Al ledt 1 

which is equivalent to the sum of the integrated mean curvatures at each of the edges, 
now weighted simply by the fraction of their associated circumcentric duals A/ that lie 
within the triangle t. A distributed mean curvature can also be defined for each triangle 
t, by spreading the integrated mean curvature evenly over the area of the triangle. 

As mentioned earlier, the edge lengths for a simplicial manifold cannot be chosen 
arbitrarily, with the triangle inequality requiring satisfaction. For the circumcentric 
dual to give a Voronoi lattice, i.e. to subdivide S into regions of points closer to a single 
vertex than to any other, the triangles must also satisfy the Delaunay condition. For 
situations where this is not true, see [7] and Appendix A of [14] . However for a surface 
embedded in E 3 , there are procedures for forming Delaunay triangulations from a given 
triangulation [15], H6J • 

The intuition used above also works best for ‘well-centered’ triangulations, where 
the circumcenter of each triangle is contained within the triangle. This is stronger than 
the Delaunay condition. A triangle which is not well-centered will have a ristricted dual 
edge which lies outside the triangle. For a consistent integrated mean curvature, this 
restricted edge must be considered to have a negative length. This leads to a reduced 
hybrid cell, again lying outside the triangle, with negative area. However, this should 
not cause much difficulty for reasonably smooth triangulations, where there is no large 
difference between the curvatures at different edges of each triangle. 

4. Higher Dimensions and Non-Euclidean Manifolds 

In this section it is shown that there is no essential difficulty in extending the 
circumcentric based distribution of the extrinsic curvature in higher dimensions, and 
more general embeddings. The generalization of a piecewise-flat simplicial manifold of 
arbitrary dimension S n , without reference to any embedding, is first described in section 
14.11 This includes higher dimensional dual lattices and hybrid cells. The distribution 
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of the extrinsic curvature for an embedding in E n+1 is given in section 14.21 requiring 
very little adjustment from the 2-dimensional case. More general situations are then 
discussed briefly in section 14.31 

4-1. Higher dimensional simplicial manifolds 

In higher dimensions, the concept of a triangle is generalized to an n-simplex. In 
Euclidean space these are n-volumes, formed by n + 1 vertices, with straight lines joining 
each pair of vertices. In one, two and three dimensions, n-simplices are given by line 
segments, triangles and tetrahedra respectively. The boundary of each n-simplex is 
formed by a collection of (n — l)-simplices, n + 1 in number. For example a triangle is 
bounded by 3 line segments, and a tetrahedron by 4 triangles. These (n — l)-simplices 
must also satisfy a generalization of the triangle inequality, to ensure that the simplex 
is well defined, and that its n-volume is non-zero. The great advantage of n-simplices 
is that their geometry in Euclidean space is entirely determined by the lengths of their 
edges. 

An n-dimensional piecewise-flat simplicial manifold S n is a manifold which can 
be decomposed into a collection of flat n-simplices, connected to each other by their 
boundary (n — l)-simplices. These are the n-dimensional equivalent of approximating 
a 1-dimensional curve by straight line segments. Since each n-simplex is determined 
entirely by its edge lengths, the internal geometry of S n is again determined completely 
by the lengths of all of its edges l. 

The Euclidean geometry interior to each n-simplex can again be developed across 
its (n — l)-simplices to any adjacent n-simplex. The deviation of the internal geometry 
of S n from Euclidean space now shows up in the deficit angles around an (n —2)-simplex. 
The circumcentric dual lattice is formed by edges connecting the circumcenters of the 
n-simplices, passing perpendicularly through the circumcenters of the (n — l)-simplices. 
These form the ‘hinges’ h on which the extrinsic curvature of S n is concentrated. The 
corresponding circumcentric dual edge to an (n — l)-simplex h is denoted A h- 

The n-dimensional hybrid cell Hh is formed by joining the vertices of h with those 
of X h , forming a complete tessellation of S n . An orthonormal basis can again be set up 
within each Hh, consisting of a unit vector along A h, a unit vector along the line joining 
the circumcenter of h to the circumcenter of one of the (n —2)-simplices bounding it, and 
continuing until an edge l is reached. This can be seen in figure |4] for the 3-dimensional 
case, also see mi for more details. This again defines a piecewise-linear Cartan frame, 
or n-bein on S n , with some degeneracy in the choice of unit vectors interior to h (in 
figure [4] there are three choices for the pair rri] and l). Due to the orthogonality of h 
and Xh, the volume of H h can easily be computed: 

H h = — h X h . (13) 

n 

The simplicity of this computation is a major advantage of these hybrid cells. 

Care should again be taken to ensure that S n gives a Delaunay triangulation, 
although this gets more difficult with increasing dimension. Recent progress has 
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(a) Two tetrahedra and their hybrid cell 


(b) Hybrid cell (A h direction expanded) 


Figure 4. Two tetrahedra, with the hinge h given by the triangle separating them, 
and the dual edge A h joining the circumcenters of the two tetrahedra ci and C 2 , (a). 
The dual edge can also be seen to pass through the circumcenter Ch of h. In (b), the 
average distance A ave across Hh perpendicular to h is displayed. The line mi joining 
the circumcenter of the hinge Ch to the center of the edge l is also shown. This line 
is perpendicular to both Ah and l. The triplet (l, mi, Ah) thus form an orthonormal 
frame within Hh- 


been made in perturbing given triangulations to form Delaunay triangulations on non- 
Euclidean manifolds however, irrespective of any embedding pn m- Again, the 
intuition in the following relies on circumcenters lying in the interior of n-simplices, 
however reasonable deviations should not cause problems. 

4-2. Simplicial hypermanifolds in E n+1 

For an n- dimensional simplicial manifold S n , embedded in a Euclidean space of one extra 
dimension E n+1 , there is a unique normal vector to each n-simplex. The arguments for 
the 2-dimensional case apply in the same way here. The angle 6 between the normal 
vectors on either side of an (n — l)-simplex h, occur only in the direction orthogonal to 
h , i.e. in the direction of the dual edge Ah- The integral of the mean curvature K from 
one side of h to the other, is given by the integral of the Kf component of the curvature, 
and equal to 6. Defining a convex region D , enclosing h and no other (n — l)-simplex, 
the integrated mean curvature over D is: 



(14) 


with the inner integral in the second part taken along each path 7 parallel to Ah- 

As in the 2 -dimensional case, this equation is unchanged in the limit as D converges 
on h , so that the integrated mean curvature along a hinge can be taken as: 


(15) 


IMC h = he , 


generalizing ([5]) to higher dimensions. Instead, however, taking the region D as the 
circumcentric hybrid cell Hh seems to give a much richer interpretation. A distributed 
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mean curvature is given by spreading IMCh evenly over H /p 

IMC h h6 nd 


K h = 


(16) 


H h i-h \ h 

n n n 

giving a more general form of (JTJ). This can again be viewed as the integrated mean 
curvature along the average distance across H h , orthogonal to h. Integrating the 
distances along a A h direction within Hh, over the (n — l)-simplex h, gives the volume 
of the hybrid cell Hh itself. The average distance d ave is then found by dividing again 
by the volume of h, giving a densitised mean curvature of: 

9 9 6 n9 

A h 


K h = 


da 


H h /h i A, 


(17) 


using (TT3|) for the volume of H h . This is equivalent to (fT6j) above, and generalizes (JHJ) , 
giving a discrete analogue of the infinitesimal extrinsic curvature. 

The integrated curvature over an n-simplex s can be found using fl9]) , (TTOl) and (TIT]) , 
with the appropriate n-dimensional substitutions, giving: 


IMC S = K h R h\s , 
h£ds 


(18) 


with the reduced hybrid cell Hh\ s representing the restriction of Hh to that part which 
lies in s. The result of (1T2]) also generalizes, from (fl3|) and (ITS]) : 

IMC S = Y]h9h^ = y ^ IMCh , (19) 

h£ds h h£ds h 

with the integrated mean curvature over an n-simplex s given by the sum of the 
integrated mean curvatures associated with each of its bounding (n — l)-simplices h, 
from (IT5j) . weighted according to the fraction of the corresponding dual edges A h which 
lie in s. 


f.3. More general submanifolds 

The definitions can be generalized for submanifolds of non-Euclidean manifolds M in 
the usual way. The simplicial manifold S n may be considered as an approximation of a 
smooth submanifold H of M, with the edge lengths in S n defined by geodesic lengths 
in H , and parallel transport of normal vectors requiring a connection on M. 

The simplicial manifold S n may also be embedded in a larger piecewise-linear 
simplicial manifold S n+m . In this case each simplex in S n lies on the boundary of 
some (n + l)-simplex in S n+m . The results of section EPl can be directly applied here, 
due to the flat space on the interior of each (n+ l)-simplex, and its development around 
a hinge of S n . 

For submanifolds with a higher codimension, for each n-simplex s there is a unique 
normal vector h s from which the extrinsic curvature can be defined. The complete 
embedding of S n will also require information about the change of span{T s S,h s } across 
a hinge h. However this is no more complicated than a 1-dimensional curve in E 3 . as 
described using the Seret-Frenet frame, see [3] for example. 
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5. Conclusion 

The circumcentric hybrid cells Hy have been shown to give a very natural volume for 
distributing the local integrated mean curvature. This gives a pointwise mean curvature 
for each hybrid cell, which completely determins the extrinsic curvature within each cell. 
These circumcentric hybrid cells have a number of clear advantages over other volumes: 
(i) they provide a complete tessellation of the simplicial manifold, (ii) a geometric 
orthonormal basis is provided by their circumcentric construction, giving a particularly 
straight-forward expression for the volume (fT3|) . and (Hi) the average distance parallel 
to the dual edge gives the precise distance for a fractional rate of change of the normal 
vectors (Eli- 

The use of circumcentric hybrid cells for the distributed mean curvature fits well 
with a broader use of the circumcentric dual lattice in Regge calculus. The circumcentric 
polytopes give a natural area for distributing the Gaussian curvature, and the other 
circumcentric hybrid cells have been shown to give consistent definitions for the Riemann 
and Ricci tensor, and the Riemann scalar curvatures [9, 10]. The circumcentric dual 
lattice arrises naturally in the Null Strut Calculus, where light rays (null struts) from 
the vertices of each spacelike tetrahedron meet along the timelike normal based at the 
circumcenter of the tetrahedron [8]. It is also consistent with the use of circumcentric 
dual lattices in Discrete Differential Geometry, and particularly in Discrete Exterior 
Calculus [19] |20j. 
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